Network analysis between neuron dysfunction and neuroimmune response based on neural single-cell transcriptome of COVID-19 patients.

Despite global vaccination efforts, COVID-19 breakthrough infections caused by variant virus continue to occur frequently, long-term sequelae of COVID-19 infection like neuronal dysfunction emerge as a noteworthy issue. Neuroimmune disorder induced by Inflammatory factor storm was considered as a possible reason, however, little was known about the functional factors affecting neuroimmune response to this virus. Here, using medial prefrontal cortex single cell data of COVID-19 patients, expression pattern analysis indicated that some immune-related pathway genes expressed specifically, including genes associated with T cell receptor, TNF signaling in microglia and Cytokine-cytokine receptor interaction and HIF-1 signaling pathway genes in astrocytes. Besides the well-known immune-related cell type microglia, we also observed immune-related factors like IL17D, TNFRSF1A and TLR4 expressed in Astrocytes. Based on the ligand-receptor relationship of immune-related factors, crosstalk landscape among cell clusters were analyzed. The findings indicated that astrocytes collaborated with microglia and affect excitatory neurons, participating in the process of immune response and neuronal dysfunction. Moreover, subset of astrocytes specific immune factors (hinged neuroimmune genes) were proved to correlate with Covid-19 infection and ventilator-associated pneumonia using multi-tissue RNA-seq and scRNA-seq data. Function characterization clarified that hinged neuroimmune genes were involved in activation of inflammation and hypoxia signaling pathways, which could lead to hyper-responses related neurological sequelae. Finally, a risk model was constructed and testified in RNA-seq and scRNA data of peripheral blood.

Network analysis between neuron dysfunction and neuroimmune response based on neural single-cell transcriptome of COVID-19 patients

Despite global vaccination efforts, COVID-19 breakthrough infections caused by variant virus continue to occur frequently, long-term sequelae of COVID-19 infection like neuronal dysfunction emerge as a noteworthy issue. Neuroimmune disorder induced by Inflammatory factor storm was considered as a possible reason, however, little was known about the functional factors affecting neuroimmune response to this virus. Here, using medial prefrontal cortex single cell data of COVID-19 patients, expression pattern analysis indicated that some immune-related pathway genes expressed specifically, including genes associated with T cell receptor, TNF signaling in microglia and Cytokine-cytokine receptor interaction and HIF-1 signaling pathway genes in astrocytes. Besides the well-known immune-related cell type microglia, we also observed immune-related factors like IL17D, TNFRSF1A and TLR4 expressed in Astrocytes. Based on the ligand-receptor relationship of immune-related factors, crosstalk landscape among cell clusters were analyzed. The findings indicated that astrocytes collaborated with microglia and affect excitatory neurons, participating in the process of immune response and neuronal dysfunction. Moreover, subset of astrocytes specific immune factors (hinged neuroimmune genes) were proved to correlate with Covid-19 infection and ventilator-associated pneumonia using multi-tissue RNA-seq and scRNA-seq data. Function characterization clarified that hinged neuroimmune genes were involved in activation of inflammation and hypoxia signaling pathways, which could lead to hyper-responses related neurological sequelae. Finally, a risk model was constructed and testified in RNA-seq and scRNA data of peripheral blood.

MutCov: A pipeline for evaluating the effect of mutations in spike protein on infectivity and antigenicity of SARS-CoV-2.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causing an outbreak of coronavirus disease 2019 (COVID-19), is a major threat to public health worldwide. Previous studies have shown that the spike protein of SARS-CoV-2 determines viral infectivity and major antigenicity. However, the spike protein has been undergoing various mutations, which bring a great challenge to the prevention and treatment of COVID-19. Here we present the MutCov, a pipeline for evaluating the effect of mutations in spike protein on infectivity and antigenicity of SARS-CoV-2 by calculating the binding free energy between spike protein and angiotensin-converting enzyme 2 (ACE2) or neutralizing monoclonal antibody (mAb). The predicted infectivity and antigenicity were highly consistent with biologically experimental results, and demonstrated that the MutCov achieved good prediction performance. In conclusion, the MutCov is of high importance for systematically evaluating the effect of novel mutations and improving the prevention and treatment of COVID-19. The source code and installation instruction of MutCov are freely available at http://jianglab.org.cn/MutCov.

DLpTCR: an ensemble deep learning framework for predicting immunogenic peptide recognized by T cell receptor.

Accurate prediction of immunogenic peptide recognized by T cell receptor (TCR) can greatly benefit vaccine development and cancer immunotherapy. However, identifying immunogenic peptides accurately is still a huge challenge. Most of the antigen peptides predicted in silico fail to elicit immune responses in vivo without considering TCR as a key factor. This inevitably causes costly and time-consuming experimental validation test for predicted antigens. Therefore, it is necessary to develop novel computational methods for precisely and effectively predicting immunogenic peptide recognized by TCR. Here, we described DLpTCR, a multimodal ensemble deep learning framework for predicting the likelihood of interaction between single/paired chain(s) of TCR and peptide presented by major histocompatibility complex molecules. To investigate the generality and robustness of the proposed model, COVID-19 data and IEDB data were constructed for independent evaluation. The DLpTCR model exhibited high predictive power with area under the curve up to 0.91 on COVID-19 data while predicting the interaction between peptide and single TCR chain. Additionally, the DLpTCR model achieved the overall accuracy of 81.03% on IEDB data while predicting the interaction between peptide and paired TCR chains. The results demonstrate that DLpTCR has the ability to learn general interaction rules and generalize to antigen peptide recognition by TCR. A user-friendly webserver is available at http://jianglab.org.cn/DLpTCR/. Additionally, a stand-alone software package that can be downloaded from https://github.com/jiangBiolab/DLpTCR.

N439K Variant in Spike Protein Alter the Infection Efficiency and Antigenicity of SARS-CoV-2 Based on Molecular Dynamics Simulation.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causing an outbreak of coronavirus disease 2019 (COVID-19), has been undergoing various mutations. The analysis of the structural and energetic effects of mutations on protein-protein interactions between the receptor binding domain (RBD) of SARS-CoV-2 and angiotensin converting enzyme 2 (ACE2) or neutralizing monoclonal antibodies will be beneficial for epidemic surveillance, diagnosis, and optimization of neutralizing agents. According to the molecular dynamics simulation, a key mutation N439K in the SARS-CoV-2 RBD region created a new salt bridge with Glu329 of hACE2, which resulted in greater electrostatic complementarity, and created a weak salt bridge with Asp442 of RBD. Furthermore, the N439K-mutated RBD bound hACE2 with a higher affinity than wild-type, which may lead to more infectious. In addition, the N439K-mutated RBD was markedly resistant to the SARS-CoV-2 neutralizing antibody REGN10987, which may lead to the failure of neutralization. The results show consistent with the previous experimental conclusion and clarify the structural mechanism under affinity changes. Our methods will offer guidance on the assessment of the infection efficiency and antigenicity effect of continuing mutations in SARS-CoV-2.

Global characterization of B cell receptor repertoire in COVID-19 patients by single-cell V(D)J sequencing.

The world is facing a pandemic of Corona Virus Disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Adaptive immune responses are essential for SARS-CoV-2 virus clearance. Although a large body of studies have been conducted to investigate the immune mechanism in COVID-19 patients, we still lack a comprehensive understanding of the BCR repertoire in patients. In this study, we used the single-cell V(D)J sequencing to characterize the BCR repertoire across convalescent COVID-19 patients. We observed that the BCR diversity was significantly reduced in disease compared with healthy controls. And BCRs tend to skew toward different V gene segments in COVID-19 and healthy controls. The CDR3 sequences of heavy chain in clonal BCRs in patients were more convergent than that in healthy controls. In addition, we discovered increased IgG and IgA isotypes in the disease, including IgG1, IgG3 and IgA1. In all clonal BCRs, IgG isotypes had the most frequent class switch recombination events and the highest somatic hypermutation rate, especially IgG3. Moreover, we found that an IgG3 cluster from different clonal groups had the same IGHV, IGHJ and CDR3 sequences (IGHV4-4-CARLANTNQFYDSSSYLNAMDVW-IGHJ6). Overall, our study provides a comprehensive characterization of the BCR repertoire in COVID-19 patients, which contributes to the understanding of the mechanism for the immune response to SARS-CoV-2 infection.

Comprehensive analysis of TCR repertoire in COVID-19 using single cell sequencing.

T-cell receptor (TCR) is crucial in T cell-mediated virus clearance. To date, TCR bias has been observed in various diseases. However, studies on the TCR repertoire of COVID-19 patients are lacking. Here, we used single-cell V(D)J sequencing to conduct comparative analyses of TCR repertoire between 12 COVID-19 patients and 6 healthy controls, as well as other virus-infected samples. We observed distinct T cell clonal expansion in COVID-19. Further analysis of VJ gene combination revealed 6 VJ pairs significantly increased, while 139 pairs significantly decreased in COVID-19 patients. When considering the VJ combination of α and ß chains at the same time, the combination with the highest frequency on COVID-19 was TRAV12-2-J27-TRBV7-9-J2-3. Besides, preferential usage of V and J gene segments was also observed in samples infected by different viruses. Our study provides novel insights on TCR in COVID-19, which contribute to our understanding of the immune response induced by SARS-CoV-2.